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Using a rigorous method of matched asymptotic expansions, I derive the equation of motion 
of a small, compact body in an external vacuum spacetime through second order in the body's 
mass (neglecting effects of internal structure). The motion is found to be geodesic in a certain 
locally defined regular geometry satisfying Einstein's equation at second order. I outline a method 
of numerically obtaining both the metric of that regular geometry and the complete second-order 
metric perturbation produced by the body. 
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Introduction. The governing equation of general rela- 
tivity, the Einstein field equation (EFE), describes how 
bodies influence spacetime curvature and move within 
the resultant curved geometry. Yet, since the seminal 
work of Einstein, Infcld, and Hoffman in 1938 study 
of this nonlinear problem of motion has largely focused on 
the post-Newtonian limit of slow motion and weak fields. 
In the strong-field regime, bodies have instead typically 
been approximated as test bodies moving in a spacetime 
that is unaffected by them. Only within the last fifteen 
years 0, H| has there arisen an analytical description of 
the gravitational backreaction: the body's perturbativc 
effect on spacetime geometry and that perturbation's ef- 
fect on the body's motion. In the case of a small mass 
m, this backreaction is called the gravitational self-force, 
and it is now well understood at linear order in m 043] ■ 
Beyond its foundational role, the self-force is also poten- 
tially of great astrophysical importance, as it describes 
the evolution of extreme-mass-ratio inspirals (EMRIs), 
in which a stellar black hole or neutron star spirals into a 
supermassive black hole. Such systems are predicted to 
be key sources for the planned gravitational wave detec- 
tor LISA , and they will afford both a unique probe of 
strong-field dynamics and a map of the spacetime near 
black holes. The self-force also provides an essential point 
of comparison with other treatments of the problem of 
motion: it complements post- Newtonian theory 0, [13] 



and fully nonlinear numerical simulations jlfjl llj ] . both 
of which are ill-suited to extreme mass ratios; and it 
fixes mass-dependent parameters in Effective One Body 
(EOB) theory [TM3. 



However, to extract orbital parameters from a wave- 
form emitted by an EMRI, one requires a theoretical 
description accurate to second order in the small mass, 
as shown by either naive or rigorous fl5| scaling argu- 
ments. Furthermore, comparisons with numerical simu- 
lations suggest that the second-order self-force would pro- 
vide a highly accurate description of intermediate-mass- 
ratio binaries and even a reasonably accurate description 
of similar-mass binaries HH, both of which should 
soon be observed by the ground-based detector Advanced 
LIGO [IH 12 1 • The second-order force would also fix 



EOB parameters quadratic in m. Although some work on 
the second-order problem has been done flil . flit , it was 
performed in an impractical gauge, with no clear means 
of calculating the force or the perturbation producinj 
it, and the basis of its approach was problematic 
0, • Detweiler has recently 21 1 examined the general 
features of the second-order problem, but he made use of 
ill-defined equations and assumed (rather than derived) 
the equation of motion and a regularized second-order 
stress-energy tensor. Hartc [22[ has derived an equation 
of motion valid at all orders, but his derivation does not 
apply to the motion of black holes. None of these studies 
has provided a definite expression for, or means of calcu- 
lating, the physical second-order metric perturbation or 
the piece of it that determines the motion. 

In this letter, I present the first complete description, 
explicitly determining both the equation of motion and 
the metric perturbation. I use the self-consistent formal- 
ism presented in Q. This formalism incorporates the 
small body's finite size, involving no infinities or regular- 
ization; it defines a worldlinc 7 that reflects the body's 
motion (even a black hole's) on any timescale; it deter- 
mines the equation of motion directly from the EFE, with 
no further axioms; and it presents the perturbativc EFEs 
in a hyperbolic form convenient for numerical implemen- 
tation. For simplicity, I take the body to be spherical 
and nonspinning, neglecting higher multipole moments. 
I work in units of G = c = 1. Greek indices range from 
to 3. Lowercase Latin indices refer to spatial coordinates. 
Further details will appear in a follow-up article . 

Self- consistent formalism. I combine two approximate 
solutions to the EFE, utilizing the method of matched 
asymptotic expansions [7I [23l - |25| . Suppose g M „(e) is an 
exact solution containing the small body on a manifold 
M. , where e is an expansion parameter that counts pow- 
ers of the body's mass. Now let r be some measure of 
distance from the body and TZ represent the spacetimc's 
lcngthscales, excluding those of the body itself. For r ~ 
TZ, well outside the body, in any global coordinates x a in 
a vacuum region fi (e.g., Boyer-Lindquist coordinates of 
the supermassive black hole in an EMRI) , I use the outer 
expansion g^ v {x a 1 e) = g flu (x a ) + h fJil/ (x a ,e) on a man- 
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ifold Me- {g^vi-M-E) defines an external background 
spacetime with no small body in it, and h^ v {x,e) — 

S,i>i e "^i T /3 ( x 'i 7( e )) describes perturbations due to the 
body, whose motion in Me is represented by 7. For 
r ~ elZ, very near the body, the metric varies rapidly, and 
there, in any coordinates (T, X a ) approximately centered 
on the body, I use the inner expansion g p!/ (T, X a /e, e) = 

gi ^(T,X a /e) + 'E ri > 1 e n H$(T,X a /e) 011 a manifold 
Mi- (gifiujMi) is the body's spacetime were it isolated, 

and H$ describes perturbations due to interactions with 
the external spacetime. The scaled coordinates X a /e 
serve to keep the body's mass and size fixed in the limit 
e — > 0, sending all other distances toward infinity; the 
use of a single rescaling factor makes the approximation 
most appropriate for compact bodies, in which the lin- 
ear dimension is comparable to the mass (in geometrized 
units). Scaling only distances, not T, is equivalent to 
assuming the body possesses no fast internal dynamics. 

In a buffer region around the body, defined by e <C 
r/lZ <C 1, either expansion may be used, and since they 
approximate the same metric, they must agree: the inner 
expansion can be expressed in unsealed coordinates and 
then expanded for r/lZ S> e, the outer expansion can be 
expanded for r/lZ <C 1, and the two results must match 
order by order in r and e. It follows that hffi ~ 
l/r n + 0(r~ n+1 ), and the 1/r™ term is determined by 
the (n — l)th multipolc moment of gip V . If r is a radial 

(n) 

coordinate centered on 7 and this form of h^J holds true, 
then both the body — in the full spacetime — and 7 — in 
the background spacetime — lie in the region surrounded 
by the buffer. If all mass dipole terms also vanish in 
this coordinate system, then the body is appropriately 
centered on 7. 

In standard perturbation theory, the linearized EFE 
would constrain 7 to be a geodesic in g^ v . To avoid this, 
I write the full, nonlinear EFE in 51 in relaxed form by 
imposing the Lorcnz gauge condition on the whole of h^v, 
rather than on each term hffl : 

g vp ^ph illJ = 0, where = - ^g^g pa h pa . (1) 

The exact vacuum EFE in il, = 0, then splits into a 
sequence of wave equations, 

E^[h^}=0, (2) 
E lu ,[hW] = 2pR lu ,[hW], (3) 

etc., where Ep U is the linear wave-operator Ep U [h] = 
n/jfu/ + 2i? Al p 1 , cr /i o - |0 , R^ypa is the Riemann tensor of the 
external background g^, and 6 2 R fJiL , is the part of the 
Ricci tensor quadratic in the metric perturbation. No 
stress-energy tensor for the body appears here, since the 
body lies outside Q. Equations ([2]) and ((3|) can be solved 
for arbitrary 7, the equation of motion of which is then 
determined by the gauge condition. That condition will 



involve 7's e-dependent acceleration o m (t, e), where r is 
proper time on 7, and I split Eq. JT]) into a sequence of 
equations for each hffl by substituting into it an expan- 
sion = E„>o e " a( " )M - 

Outer expansion. In the buffer region, I construct the 
most general possible solution to the wave equations ([5]) 
and ([3]) and gauge condition (TTJ. That local solution 
is then used to construct a global solution. I work in 
Fermi- Walker coordinates (t,x a ) centered on 7, where t 
coincides with r on 7 and x l are Cartesian coordinates 
on the spatial submanifold transverse to 7 at time t. I 
assume each h^J can be expanded for small geodesic dis- 
tance r = y/ 5ijX l xf So at first order, 

E r m h^ e \t)n L , (4) 

m>-l,l>0 

where L = i\ ■ ■ ■ it, n l = x l jr is a radial unit vector, and 
n L = n^ 1 ■ ■ ■ n 1 ^ is the symmetric and trace-free (STF) 
combination of t such unit vectors; this decomposition 
in terms of n L is equivalent to an expansion in spherical 
harmonics [26]. Substituting (Q| into ([2]), one finds that 
each term in Eq. Q must satisfy a Poisson equation, 

E/a^™^)^'^*) = ^- 2 E^t m, ' ) W» i The 

source S^™ is a linear combination of the lower-order 

terms hr^™, <m ' 1 it involves derivatives, a' 1 , and Rp V pa- 
The general solution to this Poisson equation consists 
of a homogeneous solution, comprising the single mode 
^(i,m,e H ) w T [ieTe £ — m f or m > and £h = —rn — 1 for 
m < 0, plus a particular inhomogeneous solution with 
modes P directly proportional to S^™' 1 \ Given 

how S}j£ is constructed, clearly each h^™'^ is a lin- 
ear combination of the lower-order modes h Ti - m ' t ' H '- 
therefore, the general solution [of the form Q] to Eq. @ 
is wholly determined by the functions hh™^ (t). The 

first of these, hfy 1,0 \ is fixed to be 2mS fiu , where m is 
the ADM mass of the inner background gi pv [fjj . All 

the others, h^',T~°' 1 are undetermined at this stage. 
It will thus prove convenient to split the general solution 
into two pieces: h$ = h$ S + h$ R . I define the regular 
field hj}} R to comprise all terms involving the undeter- 
mined functions /if, 1 '7 i ~°'^ ff ' ) . It reads 

^ )R = ^°' C0 (t)+rfc& M) (*)n < + O(r a ). (5) 

The singular field h^i then comprises all the other 
terms; it is the particular solution obtained by setting 

C^f ' eH) = 0- » reads 

+ rJ2 h%Y\t)n L +0{r 2 ), (6) 

£=0,2 
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where the functions h^J^' e \t) arc linear in m. I stress 
that this split is merely a convenient grouping of terms in 
the general solution in the buffer region, with no impact 
on the results. When so defined, h^} S and h$ R are each 
solutions to the first-order field equations in f2. h^J can 
be interpreted as the body's bound field; it is determined 
solely by the fact that a compact body lies in the region 

( 1)R 

surrounded by the buffer. h^i is a homogeneous (local) 
solution to the wave equation even at r = 0, propagating 
independently of the body; as a homogeneous solution, it 
can be determined only by global boundary conditions. 
Imposing the gauge condition determines the isotropic, 
5^ form of the m/r term given above, as well as deter- 
mining that the body behaves approximately as a test 
particle, with constant mass (i.e., dtm = 0) and approx- 
imately geodesic motion (i.e., a^ ^ = 0). 

One can prove with distributional [j| or Green's func- 
tion Q methods that h$, because of its m/r term, is 
sourced by the stress-energy tensor of a point particle, 
T^ v = ^mu^u" Si{x "jll {T)) dr, where z»(t) is the 

parametrization of 7, = ^f— , and g is the determi- 
nant of g^y. Therefore, at distances in the buffer region 
or greater, the body appears as a point mass. With this 
determined, hfpu can be found globally by solving the 

wave equation E^^h^] = —16ttT„ v , where an overbar 
indicates trace-reversal. Using retarded boundary condi- 
tions, for example, the global solution will then fix the lo- 



cally undetermined field h$ R . Numerous methods have 

been used to accomplish this [4). In the case that h$ 
contains no contribution from incoming waves at infin- 
ity, doing so determines Q that at least through order r, 



is the Dctweilcr- Whiting regular field [27 



The second-order solution proceeds almost identically. 
I assume an expansion 

m>-2,l>0 



In r 

m>0,£>0 



m , (2,m,ln,£) 



m 1 



(7) 



where the logarithmic terms arise from the correction 
~ mln(r/2m — I) to the light cones in the buffer re- 
gion (where r ^> m) @. Substituting this into Eq. 
together with an expansion of h$ up to order r 2 , and 
finding the general solution at each order again allows 
the convenient split /ijfj = hffi R + hffi S . However, 



here I define h 



(2)R 



to comprise not just all terms in- 

but 



volving the undetermined functions h^™^ ' e "\t) 



( 1)R 

also all terms quadratic 111 tlui/ This guarantees that 
(2)R t? ru(2)Ri _ 26 2 R llv [h^ K ], such that 



h\tr satisfies E^Jh^ 



9^ 



Am 



,(2)R 



satisfies the vacuum EFE through 



order e . Explicitly, 



"'HV 



h^M+rh^HW + Oir 2 )-, (8) 



TABLE I. Components of <5m M „ (t) in terms of the first-order 
regular field h$ n (t,r = 0) = h^) m (i). 



-rah 



(1,0,0) 



- ±m<5 1J 7i (1 '°' 0) 

tt 6 ij 

(1,0,0) 



5m u 

8m ta — — |mft.__ 
Sm ab = S ab (imhg' ' 03 + ^m^hf^ 0) 



3 {ab) 



terms quadratic in hfy R would appear at order r 2 . The 

(2)S 

singular field h^J , comprising all other terms in the gen- 
eral solution, then reads 



hi 



2)S 



1=1 



e=a,2 

1=1 £=0,2,3,4,5 

+ lnr[/#°>>°> + rhft^n*] + 0(r 2 ), (9) 

where 5m M „ is a mass-like tensor defined on 7. The ex- 
plicit terms in this expansion can be found in Ref. @ 
through order r° when a M = 0, and they will be writ- 
ten out in full in the follow-up article [2(J; the functions 
(t) arc linear combinations of m 2 , mh^ R ', and 
The body's dipoles (i.e., those of gi^ v ) would 

,(2)S 



h (2,m,l) 
n HvL 

5m,,,,. 



contribute to h\i , but I set the spin to zero for sim- 
plicity and the mass dipole to zero to ensure that 7 ac- 
curately represents the body's motion. Note that my 

(2)R (2)S 

definition of h^i means that K^l satisfies not Eq. j3J), 
but E^h^ 8 ] = 25 2 R^[hW] - 25 2 R lf „[hW R ]. 

As at first order, the gauge condition determines the 
form of 5m ^ , as given in Table HI along with the first- 
order acceleration. 



(*) 



,(1,1,1) 



(1,0,0) 



(10) 



The evaluation at = is to be performed only at time 
t, leaving the past history of 7 unchanged; this follows 
from the presumed expansion of a^, and it prevents a 
need for order-reduction 0, 0] • One can show Eq. (fT0|) 

(1)R 

is equivalent to the geodesic equation in g^ + th^i at 
order e 

Using the same Green's-function method @ as at first 
order, one can straightforwardly prove that the terms 
involving 5m^ v have a point source, the trace-reversal of 
which is given by the effective stress-energy tensor 



" ///,' 



dr. 



-9 



(11) 



The global solution to the second-order EFE sourced by 
the body can then be obtained numerically via a puncture 
scheme [28[: Outside a tube T around the body, one can 
solve Eq. directly; inside the tube, one can use an 
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approximation hffi S to h$ a , given by Eq. (J9j) without 



,(2)S 



the "0(r 2 )" term, and a regular field hffv = h)£v — 
which satisfies 



,(2) I(2)S 



E. u [¥ 



2)Ri 



16 7 rT2 y ) +2 ( 5 2 J R^[/ l (1 )]-^[/ 1 ( 2 ) s ]. (12) 



All divergent terms on the right-hand side cancel, leav- 
ing fijjj, to solve a wave equation with a regular source. 
At r, the analytical expression for h\a, can be added to 
the numerical solution to Eq. (fT2"|) to change variables to 
the full field. Since hffi R will agree with hffi R through 
order r, this procedure will also determine hpu° (t) 
and hfyl' (t), which, as we shall find. This puncture 
scheme can be implemented immediately after transform- 
ing Eq. ([5]) to a desired coordinate system. 

Matching to an inner expansion. One could proceed 
with the same method to find h^} together with e^ 2 ^. I 
instead take a more efficient route by determining a*- 2 ^ 
from additional information about the inner expansion. 
I take the small body to be a Schwarzschild black hole. 
Since the inner expansion affects the outer one solely 
through the body's multipolc moments, this amounts to 
neglecting non-monopole moments, as was done in the 
preceding section. Beyond effects of those moments, the 
equation of motion I derive will hold for any compact 
body. I also specify the perturbations to be produced 
by tidal fields. These fields are of quadrupole order and 
higher, and while they produce mass and spin perturba- 
tions via tidal heating and torquing, dimensional anal- 
ysis shows those perturbations do not contribute at the 
orders in e of interest. With these specifications, I fol- 
low the procedure in Ref. 0, insisting that in a suitable 
mass-centered coordinate system, this inner expansion 
must equal the outer expansion in Fermi- Walker coordi- 
nates when expanded in the buffer region, up to a unique 
gauge transformation that excludes spatial translations 
at 7; this ensures the desired relationship between 7 and 
the mass-centered inner expansion. Here 'mass-centered 
coordinate system' means one in which the mass dipole 
°f 9i {iv vanishes along with all even-parity dipole per- 
turbations that behave as a mass dipole, scaling as 1/r 2 
in the buffer region. Other even-parity dipole perturba- 
tions are also set to zero in order to leave no residual 
gauge freedom. Because the Fermi- Walker coordinates 
are e-dependent, each term in the outer expansion de- 
pends on a M . Hence, prior to matching the metrics, I 

substitute = J2n>o e " a fn) into 9u» + L>i e n ^l" ■ 

A tidally perturbed Schwarzschild metric is given in 
Ref. [2^] in advanced Eddington-Finkelstein coordinates 
and a light-cone gauge. The tidal fields are represented 
by STF tensorial functions of time: £ij and Bij for 
electric- and magnetic-type quadrupole tides and £ijk 
and Bijk for analogous octupole tides; for the powers of 
r of interest, hexadecapole and higher tidal fiels would 
appear only at order e 3 in the outer expansion. Trans- 



forming to a suitable Fermi-like coordinate system and 
expanding the result to order e 2 in the buffer region yields 

Eu - ~fu + Mt£n^ + r 3 £ H^n L (13) 

(14) 



=2,3 



E ta = r^tfn L + r^Htt ) n L 



1=1 



(=1 



gab = f ab +fn ab+ r*Y: H^ + r^Y: H^h L . (15) 

1=0 e=o 



Here /„ = 1- ^-f 2 ^, f ab = (l+^ + ^)<U, and 
f= ! ~ 



describe the Schwarzschild metric in harmonic 

coordinates. The coefficients H^ L are functions of em/r 
forming linear combinations of £ij(t) and Bijit), and 

are l mear combinations of £ i j{t) 1 Bij(t), £ijk(t), 
and Bijk(t), where an overdot indicates a time-derivative. 



For example, H^ff 



-(1 



5 em 



+ 



4e"W 



3r 1 -iP~>^ l r 

If explicit appearances of \i are set to zero, the met- 
ric of Eqs. ([m f - pl)]) reduces to that of a vaccum space- 
time in Fermi coordinates centered on a geodesic, with 



c _ c(0) 

c _ c(0) 
°abc — O abc 



eS£ ab + 0(e 2 ), B ab = B<°> + eSB ab + 0(e 2 ), 



0(e), and B a 



be 

(0) H (0) 



B 



(0) 
abc 



O(e), where 



the zeroth-order fields £ ab , B ab , etc., are components 
of Rfiupa and its first derivative evaluated at r = 0. 



(n)R 



There is no manifest appearance of the fields h 
Eqs. (|13[) - (fT5")) ; they are incorporated into the correc- 
tions 8£ a b, 5B abl etc. More significantly, there is no term 
corresponding to an acceleration; any such term would 
induce a mass-dipolc-like term and would vanish in mass- 
centered coordinates. 

To match the expansions at orders e and e 2 , I seek a 



unique transformation x^ 



e 2 £< 2 ^ that 



brings the outer expansion into the form of Eqs. (|13|l - 
(flS). Decomposing and £ (2) '' into irreducible STF 
pieces, one readily finds a unique transformation. At 
each order in e, this transformation can be thought of 
as putting the outer expansion into Fermi coordinates in 
SV+e^uy The order-e transformation is given in Ref. 
up to order-?' terms in the metric; matching the metrics 
exhausts all freedom in that transformation and uniquely 
determines the standard result (fTU| for o^ 1 ^. Matching 
order-er 2 terms in the metric fixes S£ij and SBij in terms 
of hjtv- The order-e 2 transformation, when carried to 
order-r terms in the metric, likewise uniquely determines 



,(2) 



h 
2 



(2,1,1) 
tti 



,(2,0,0)1 



,(1,0,0) Jl) 



1 



(1,0,0), (1,0,0) 



11 



(16) 



where, as in Eq. (|T0|) . the evaluation at = occurs 
only at time t. Summing ea^-, and e 2 a^ 2 w one finds, up 



5 



to 0(e 3 ) errors, 



(17) 



where h^ v 



Am 



,(2)R 



This is the geodesic equa- 
tion in the locally defined regular metric <7 M „ + h^ v , up to 
terms cubic in h^ v . g^ v + h^ v may trivially be extended 
to a C n (local) vacuum solution to the EFE, through 
order e 2 , by finding the solution (UJ) through order r n . 

Equation (1 1 7|) agrees with the form of Harte's equation 
of motion [22[ but represents a major advance: it has 
been shown to apply to black holes, and it comes along 
with a concrete means of calculating both et^ 2 ^ and h^o ■ 
Discussion. I have shown that through second order in 
its mass, a small body moves on a geodesic of a certain 
locally defined regular metric. I have also derived results, 
given by (|12[) . and (fTrJ)) . that (together with the first- 
order equations) may be used to simultaneously evolve 
the body's position and find the perturbation due to it, 
thereby solving the EFE through second order. Although 
these results were derived only for a nonrotating black 
hole, they should hold for any spherical, compact body 
with slow internal dynamics. For nonspherical bodies, 
they will be modified by higher multipole moments: h^l 
and h)J will be straightforwardly altered by the body's 
spin 0], and a (2)/i will include well-known [H El ED] 
couplings of the moments to the external curvature. 

I wish to thank Leor Barack for many helpful discus- 
sions and suggested improvements to this manuscript. 
This work was supported by the Natural Sciences and 
Engineering Research Council of Canada. 
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